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The present status of lattice calculations of fn , /b and some mass splittings are discussed. When one includes 
the uncertainties due to discretization errors, the results do not yet have a sufficient precision to be relevant to 
phenomenological applications. There are, however, good prospects of cutting down the uncertainties by a factor 
of 2 or more soon. 



1. INTRODUCTION 

B-physics plays an important role in the 
experimental determination of the Cabibbo- 
Kobayashi-Maskawa matrix and the understand- 
ing of CP-violation. Until CP-violation can be 
observed directly in B-meson decays, hadronic 
matrix elements are needed in combination with 
the experimental results on B — B, K — K mixing 
to obtain restrictions[l] on the famous unitarity 
triangle and thus on the CP-violating phase 6. 
Only lattice QCD allows us to compute these ma- 
trix elements without model assumptions. 

In this review, we do not discuss the relevant 
B — B matrix element but concentrate on the 
leptonic decay constants /g and fjj , as well as 
beauty spectroscopy. These quantities are easier 
to compute. It is important to understand them 
before one performs a full study oi B — B mixing. 

This talk is an update of ref. [2] , which we rec- 
ommend as an introduction for the non-expert. 
We restrict our attention to the quenched approx- 
imation and refer the reader to ref. [2] for the little 
that can be said about full QCD. 

Compared to last year's conference, where the 
subject was reviewed by C. Bernard [3], there has 
not been a rapid development. However, at least 
in the static approximation, it is well understood 
how to reduce errors to around the 15% level due 
to the work of the FNAL group [4] and Draper 
and McNeile [5] as well as the variance reduction 
to be discussed in section 5. 



2. fo 

Leptonic decays, as e.g. D~ —f l-i>i, are given 
in terms of a simple vacuum to pseudoscalar ma- 
trix element 



{0\Ao{0)\P) = fp^ymp/2. 



(1) 



of the relevant axial vector current. In princi- 
ple, the leptonic decay constant fp can easily be 
obtained from a Monte Carlo estimate of the cor- 
relator {Ao{x)aI{0)). 

Nevertheless, there are two non-trivial prob- 
lems. 1) The axial vector current acquires a 
renormalization in the lattice regularization and 
2) for current values of the lattice spacing o, the 
propagation of the charm quark is distorted be- 
cause its mass is not small enough compared to 
the inverse lattice spacing. The second problem 
becomes a true obstacle for /g . 

2.1. Renormalization 

The relation between the bare lattice current 
(in the Wilson formulation) and the renormalized 
current A^'^ (a;) is given by 

Af/{x) = ZA{glKf,Kf.) qf{x)'r^'rsqf'{x), (2) 

with Kf the hopping parameter of quark flavor 
/ and Qq the bare coupling. In non-relativistic 
normalization and with tadpole improvement [7, 
8] it is natural to use 



Al'f'ix) = ZA(9n,Kf,Kf-) J— — 
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7^ qfixh^,^5qf'ix) . (3) 



2Kf- 8Kc 



Kc denotes the critical value of the hopping pa- 
rameter. These two definitions are equivalent 
when one imposes a non-perturbative normal- 
ization condition for the currents, but they are 
quite different when Za and Za are taken from 
1-loop perturbation theory neglecting terms of or- 
der aruf = \/2Kf — l/2Kc. In the numerical 
estimates, we will use Za and Za expanded in 
terms of g^ = Zg^j (TrUa), which should be a 
well-behaved perturbative series [8]. 
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Figure 1. The lattice spacing dependence of iVf^/vS 
for three different fermionic actions: Wilson (squares) SW 
(triangles) and staggered action (filled squares). Mp as 
calculated by several groups [2] had to be corrected for 
finite volume effects in some cases. S is taken from ref. [12] 
as throughout this review. Note that more recent results 
[13] at S ~ 0.16 for the SW-action are lower than the 
value shown here. 



2.2. The continuum limit 

We use capital letters for quantities in lattice 
units, e.g. Fp = afp. Since QCD predicts only 
dimensionless quantities such as mass ratios, we 
consider the continuum limit of Fp/Ft^. The 
dominant lattice artefacts should originate from 
the lowest-dimensional operator which appears in 
the effective action of lattice QCD but not in the 
continuum QCD action. For Wilson fermions, 
this is a dimension-five operator. Hence, we ex- 
pect the continuum limit to be approached at a 
rate proportional to o - up to unimportant loga- 
rithmic corrections: 



If instead we take the perturbative values for Za 
and Za and use Mp to set the scale, there are 
additional perturbative corrections 



Fp/Mp = fp/rrip + dg"" + C'Mp 



(5) 



Fp/F^ = fp/U + CF^ 



(4) 



d and higher-order perturbative terms vanish only 
when the currents are normalized through a non- 
perturbative normalization condition, e.g. the 
axial Ward identity [9]. In this case, one still 
has to remember that C" depends on which nor- 
malization condition has been chosen. In fact, 
for the case of the vector current this 0{a) ef- 
fect is known to be of the order of 15% still at 
(3 = 6.4 [10]. Furthermore, the coefficients C and 
C" depend significantly on whether one chooses 
eq. (2) or eq. (3), the difference being of order 
exp(o(my -\-mft)/2). We want to point out, how- 
ever, that in addition to the specific 0{amf) lat- 
tice artefacts of heavy mesons, there are sizeable 
generic order 0{a) effects when one uses the stan- 
dard Wilson action. To show this, we consider the 
ratio of the p-mass to the string tension in fig. 1. 
At /3 = 5.7 (VS = a-s/a ~ 0.4) there are ~ 30% 
lattice artefacts. 

The above general discussion of lattice artefacts 
serves to underline three points: 

• Perturbative corrections are potentially 
dangerous because they come as powers of 
g^ . Due to their logarithmic variation with 
the lattice spacing they can hardly be de- 
tected numerically. For the case discussed 
here, they can in principle be avoided by 
choosing eq. (4). 

• There are significant generic lattice arte- 
facts that cannot be eliminated by simply 
changing the normalization of fields. In 
order to reduce them, one needs to use 
an improved action and improved currents, 
the first candidate being the systematic 
0{a) improvement of Sheikholeslami and 
Wohlert (SW) [11]. For that action - with 
tree-level coefficient - there is indirect evi- 
dence of improvement compared to the orig- 
inal Wilson action [14]. A direct test of im- 
provement could be obtained through pre- 
cise computations of Mp/yYi in the range 
(3 = 5.7-6.0. 



• Whichever action and normalization one 
chooses, the simulation results need to be 
extrapolated to the continuum using a form 
like eq. (4). Only in this way can one 
roughly account for the uncertainty of the 
continuum results due to 0(o) lattice arti- 
facts (that are possibly hidden in the sta- 
tistical errors). 

2.3. The present status 

An extrapolation of fjj to the continuum limit 
was reported in ref. [15]. It has been repeated 
using also the results of other groups [17,19] in 
ref. [2] . Because of the relatively large statistical 
errors of F^^ in these simulations, it is convenient 
to extrapolate separately i^^/yS and Fo/yS to 
the continuum and then take their ratio. In this 
way, the 0[g'^) errors should roughly cancel and 
the results of many groups for F^^ can be taken to 
obtain the first ratio. ^ 
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Figure 2. Continuum extrapolation of i^7r/vS (squares; 
data from many groups [2] ) , Fj) / v S in relativistic normal- 
ization (filled squares [15,19,17]) and in non-relativistic 
normalization (triangles) using 7712 as described in the 
text. For comparison we show also i^^/vS with the SW- 
action (filled triangles [16,17]). 



As seen in fig. 2, the extrapolation of i^^/yS 
is rather unproblematic. Fo/yS does, however. 

Of course, the results for the decay constants need to 
be extrapolated in the mass of the light quark and an 
interpolation in the heavy quark mass needs to be done 
in order to keep e.g. iVfp/vS = mj^ j \/a at each value of 
the lattice spacing. 



depend significantly on the normalization of the 
fields. Compared to the relativistic normaliza- 
tion the data shows much stronger lattice spac- 
ing dependence in the non-relativistic normaliza- 
tion, see fig. 4 of ref. [15]. The reinterpreta- 
tion of the Wilson action as an effective non- 
relativistic action for large ma suggests [7] to 
take as the mass of the meson rap the kinetic 
mass m2 instead of the pole mass m\ in the non- 
relativistic expansion of the energy of the meson 
£^(p) = mi +p^/(2m2) + .... Since the simulations 
mostly do not compute m^-, one needs to express 
m2 in terms of m\ (and possibly the bare quark 
mass [19]) approximately. This introduces an am- 
biguity in the procedure, but generally the depen- 
dence on the lattice spacing is reduced again. The 
open triangles in fig. 2 illustrate the weak lattice 
spacing dependence for the case where one uses 
the relation between ra^a, and m\a that is ob- 
tained for a free Wilson quark. We emphasize 
that these modifications are irrelevant in the con- 
tinuum limit and one should have results over a 
large enough range of o such that after the ex- 
trapolation the details do not matter any longer. 

Inserting the experimental value of Z^, the 
above extrapolations yield /d = 164(20) MeV 
(relativistic normalization) and /d = 188(22) 
MeV (non-relativistic normalization). A safe es- 
timate is therefore Jd = 176(34) MeV 

At this conference, two new computations of 
fjj were reported. We show them in fig. 3 to- 
gether with previously published values. The pre- 
liminary results from the MILC collaboration [20] 
have very small errors, given the statistical en- 
semble. If the errors of the final analysis are of 
this order, one can indeed extrapolate directly 
Fd/Ft and avoid the uncertainties due to the 
perturbative renormalization. Unfortunately the 
precision oi Fd/F^^ and the range in lattice spac- 
ings where this quantity exists is not sufficient yet 
to establish its lattice spacing dependence with 
the SW-action (cf. fig. 3). 

It is quite obvious that /g cannot be computed 
in this way. A possible approach is given by a 
non-relativistic treatment of the b-quark [6,7]. In 
these approaches it is necessary to estimate the 
systematic errors due to uncertainties in the co- 
efficients in the action and due to the truncation 
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Figure 3. Preliminary results for Fd / F„ with the Wil- 
son action (non-relativistic normalization) from ref. [20] 
(squares), together with the numbers from ref. [21] (dia- 
mond) and [17] (triangles). Full symbols are for the SW- 
action, with the inverted triangle point from ref. [16]. 



very difFcult to obtain significant correlators be- 
yond a distance of 1/2 fm. In this situation, it is 
essential to use smearing to suppress the excited 
states. 

3.1. Computing Ground State Properties 

Concerning this problem, a true breakthrough 
has been achieved by the FNAL group [4] and 
subsequently by Draper and McNeile [5], who 
showed that the variational approach [23] can be 
applied successfully to this case. I briefly out- 
line the main principle but not the details of the 
analysis performed by the two groups. 

One constructs a matrix of correlators 



c"M = E(4''(»'o,^i(4'V(o)> 



(8) 



of the action. Applications are presented in [30]. 
We here discuss the following approach instead. 
One computes the combination 



/ = fpVmp[(Xsi'mp)/as{mB 



l2//3o 



(6) 



for a range of masses, say up to rUp ~ 3 GeV. 
In addition, Eichten's static approximation [22] 
allows one to determine 



r 



lim fp^/mp . 

mp — ^oo 



(7) 



In both parts of the calculation, one can take the 
continuum limit (after renormalizing the current 
in the static approximation). In the continuum, 
one then matches the two parts through one or 
two phenomenologically determined l/rrXp correc- 
tion terms. ^ 

3. /s IN THE STATIC APPROXIMA- 
TION 

A calculation of ^**''* is essential for the suc- 
cess of the above approach. Unfortunately, the 
static approximation suffers from a bad signal- 
to-noise ratio of the correlation functions: it is 



Such a matching cannot be justified at a finite value 
of the lattice spacing [19,20]. The reason is that in the 
regime aruf >> 1 one is essentially in the non-relativistic 
limit, where one cannot take the continuum limit but must 
add higher-order operators in order to cancel cutoff'eff'ects. 
Therefore, one may not interchange the interpolation and 
the continuum limit . 



I, J = 0, ..., M, of axial vector currents with dif- 
ferent (Coulomb gauge) wave functions 

a''J\xo,x) = ^qstatixo,x + y)joj5^^^\y) (9) 
y 

qf{xo,x-y) . 



Its spectral decomposition reads 



C"{xo) = J2<i<y^M-EcXo) 



(10) 



a=0 



With / = corresponding to the local axial vec- 
tor current, one has #'*"* = ^/2Z'*''*v^, where 
^stat jg ^j^g renormalization of the axial current in 
the static approximation [24]. Diagonalizing the 
matrix C-'^'^{to)C{t)C-'^'^{to) for t > to, yields 
estimates i5^ and Ea, which agree with the exact 
overlaps v^^ and energies Ea up to corrections of 
order 0(exp[-(£;M+i - Eo,)t]) [25]. So, one ob- 
tains a good estimate of the gap AE = Ei — Eq- 

Due to their construction, the trial wave func- 
tions $(^' of ref. [4] have a large overlap to the 
lowest M states. Ref. [5] constructs a complete 
basis of functions $. Both groups finally do not 
take j'**"* from Vq directly but use Vq to construct 
smeared-smeared and smeared-local correlation 
functions to be analyzed in the standard way. 

In fig. 4, we show that the results [4,5] are in 
good agreement. Previous values such as [26] ap- 
pear rather high (see fig. 20 of ref. [4]). Fig. 5 
demonstrates explicitly (using the estimate for 
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Figure 4. f'*""/^^''^ with the Wilson action (squares 
[4], triangle[5]) and the SW action (filled circles [16], filled 
triangle [18]). Note that for the Wilson action, we use the 
mean-field improved 1-loop values for Z^^°'^ [24]. For the 
SW action this calculation has not been done yet and the 
1-loop results in g^ are used. This is the main reason for 
the diff'erence in the results. 



the gap A£^ from ref. [4], and the tables of 
[26]) that the difference between the results of 
ref. [5] and [26] is due to an unresolved correction 
of order 0(exp(— A£^ ^min)) in the latter refer- 
ence, (the gap A£^ not being known, the results 
were considered as asymptotic within errors at 
exp(— A£^ tmin) — 0.3). Notice in particular that 
the 0(exp(— A£^ imin)) correction becomes larger 
with decreasing quark mass. Therefore, fs, / fs^ 
is strongly affected. In the upper part of the fig- 
ure, one of the results of [4] receives confirmation 
from [31]. We conclude that the difference visible 
in fig. 20 of ref. [4] is quantitatively understood 
to be due to an excited state contamination. 

Ref. [4] investigates carefully also other sources 
of systematic errors like finite-size effects, and the 
extrapolation to the physical light quark masses. 

3.2. The Continuum Limit 

Finally, F'*"* /{ll^mfl^ has been extrapolated 
[4] to the continuum limit, with Am the IP-IS 
charmonium splitting [27]. Am is not known for 
(3 = 6.3. Ref. [4] approximates it by 1-loop evolu- 
tion from /3 = 6.1 to /3 = 6.3. This does, however, 
give an arbitrary value for the o-effects of the ra- 
tio F^*'^*/{Am)^/'^ (e.g: the same procedure ap- 
plied also to the numerator simply says that the 
ratio is independent of the lattice spacing). Fur- 
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Figure 5. Dependence of the bare f"*'"' on the minimum 
time distance tmin in the fit to the local-smeared corre- 
lation function. Filled circles correspond to gauge invari- 
ant Gaussian smearing functions, while filled triangles are 
from exponentials [26]. At exp( — A.E tmin) = 0, we show 
the results of the variational calculations for comparison 
[4,5]. At /3 = 5.7 no fit to the local— smeared correlation 
function was done, but only one value t = tmin was used 
[31]. The smearing function was not optimized in that 
case. 



thermore, the statistical errors of Am have not 
been included. We therefore performed the ex- 
trapolation of p^tat ^-^3/4: ysijjg also the data of 
[5]. The lattice spacing dependence of this ratio is 
weak. ■* Conservatively, using only the last four 
points to extrapolate, gives the continuum ratio 
fstatj^zli ^ 1.85(37) with an additional (esti- 
mated) 7% uncertainty due to the renormaliza- 
tion [24]. 

In fig. 4, also the available values with the SW 
action are displayed. Since these calculations do 
not use the variational method, we suggest that 
they should be analyzed as shown in fig. 5, to 
make sure that 0(exp (— A£^ tmin)) errors are un- 
der control. Nevertheless, it is exciting to see that 
one will soon be able to obtain the continuum 



^Note that figs. 1 and 4 mean that in contrast F'^'^^/M^^^ 
decreases strongly with decreasing a, explaining the very 
high numbers that were obtained originally [28,29]. 



limit for this difficult quantity with two different 
actions. The main uncertainty that will remain 
is the 2-loop uncertainty in the renormalization. 

4. fB 

One may now combine f^tat ^^3/4 ^j^jj mpj \/a 
or /tt/v^ to obtain j^^°-^ in physical units. To- 
gether with f(mp) in the range 1.2 GeV < mp < 
3 GeV obtained as in sect. 2.2, one finally deter- 
mines /b = 180(48) MeV through interpolation. 
First results using the actions [6,7] are in agree- 
ment with this value [30]. 

A significant improvement of this result should 
be possible. Reducing the lattice spacing will al- 
low for a more precise extrapolation of f{rnp) 
and a wider range of mp. At the same time, it 
is necessary to improve the statistical accuracy, 
especially in the static approximation. In the fol- 
lowing section, we demonstrate that a significant 
factor can be obtained without additional com- 
putational effort. 

5. VARIANCE REDUCTION 

Consider for simplicity the correlation function 
of two pseudoscalar densities. After integration 
over the Grassmann variables, it can be written 

Oi{xo) = '^{TrSf{xo,x;0,y)SJ,{xo,x;0,y)), 

X 

where the average <> is over the gauge fields 
with the appropriate weight including possibly 
the fermion determinant and Sf is the quark 
propagator of flavor /. The variance of this cor- 
relation function could be decreased by averag- 
ing over y. Straightforwardly this is not pos- 
sible, since it requires quark propagators to be 
calculated from each point y. Instead, with just 
the effort necessary to compute S from one point 
y, we can calculate the combination Sf(xo,x) = 
S« 5y (so, k; 0, y)cr^ where cr^ is a random field of 
Ising variables. With this building block, a sec- 
ond observable (N = ^^ 1) 

O2ixo) = ■j;^'^{Tr Sf{xo,x)SJ,{xo,x)) 

X 

can be constructed, which on average is exactly 
equal to Oi. Here it is understood that we aver- 



age also over the Ising field. The variance of O2, 
1^(02), involves in addition to the terms that are 
present in V(Oi) 4-point functions that are partly 
summed over their arguments. If these 4-point 
functions decay fast enough as the arguments are 
separated, V{02) = 0{1/N) x V{Oi). The pref- 
actor may however also be large, such that on 
a lattice of (1.5 fm)^ one has V{02) > V{Oi) 
because N is not large enough. This prefactor 
originates from the sum over short distances in 
the 4-point functions. It can be decreased by not 
summing over every point y, but over "well sepa- 
rated" points only. 

The essential idea to reduce the variance of 2- 
point functions like this was given in ref. [32]. 
There, and in a recent investigation [33], it was 
concluded that the method does not improve the 
variance in the cases of practical interest. How- 
ever, we do not expect a large prefactor when at 
least one quark flavor is heavy. 

This variance reduction can be applied for any 
action of the quarks and for any type of smearing. 

The idea has been tested both in the static ap- 
proximation and for a heavy quark around the 
charm-mass at /3 = 5.7 on a 12^ x 24 lattice 
and with smearing [31]. As suggested by the 
above argument, one finds a shallow minimum 
of 1^(02) as a function of the separation between 
the points y. The minimum occurs around a sep- 
aration of 0.3 fm [31]. It is considerably lower 
than V{Oi): V{0i)/V{02) = 3.4 - 4.5 in the 
relevant time interval in the static approximation 
and V{0i)/V{02) = 2.6 - 4.0 for a charm-light 
correlator. One clearly expects that these ratios 
will grow proportionally to the space-like volume, 
when the separation of points is kept fixed. It was 
also checked that the gain translates e.g. into a 
factor 1/2.5 in the error on #'*"* on the 12^ x 24 
lattice. 

Such factors should not be missed in future 
heavy-light and heavy-heavy computations. 

6. BEAUTY/CHARM SPECTROSCOPY 

The spectroscopy of mesons and baryons with 
b- or c-content, is of twofold interest. On the one 
hand, there are still channels where lattice gauge 
theories can make predictions; on the other hand. 



one can check the importance of systematic er- 
rors like the quenched approximation against ex- 
perimental numbers and one can test the size of 
1/m corrections to the heavy quark limit in these 
quantities. Two splittings have been studied sys- 
tematically. 




Figure 6. Continuum extrapolation of (M^, — 

Mb^)/v^- Data are from [4] (squares) and [13] (triangle). 



The splitting Mb, — Mb^ was computed in the 
static approximation for a range of lattice spac- 
ings by ref. [4]. My remark concerning the con- 
tinuum limit extrapolation in sect. 3.2 applies 
here as well. I have done the extrapolation of 
{Mb, - MbJ/V^, see fig. 6. I have included the 
result of [13] in this case, since it is well known 
that the excited state contributions are much less 
relevant for the extraction of the masses than for 
the determination of the decay constants. Indeed, 
a preliminary value [5] obtained with the varia- 
tional technique is in agreement with ref. [13]. 

In the continuum limit one obtains {niB, — 
rriBj/Vo^ = 0.19(2) + 0(v^/ms). This is to be 
compared to the "experimental" result {niB, — 
-ttibJ/V^ = 89(4) Mey/(420 - 440) Me^ = 
0.20(1) — 0.21(1), where the real uncertainty in 
assigning a value to the string tension is hard to 
quantify. Unless there is an accidental cancella- 
tion of the two effects, the influence of dynami- 
cal fermions and the 0(l/m;,) terms are not very 
important. This is in agreement with the experi- 
mental fact that the splitting mj}, — mj}^ is only 
~ 10% higher. 

Maj — Mb has been computed with the Wilson 



action for a range of lattice spacings with vary- 
ing masses of the heavy quark [34]. The lattice 
spacing effects of (Ma^ —Mb)/Mp are small. One 
can extrapolate to the continuum and to the mass 
of the b-quark [34]. In the static approximation 
there are estimates at /3 = 5.74 [26] and at /3 = 6.0 
[35] which, however, are obtained from only mod- 
erately long plateaux. Stella has presented results 
from UKQCD obtained at /3 = 6.2 with the SW 
action [36]. They are in nice agreement with [34], 
thus indirectly confirming the smallness of lattice 
spacing effects in (Ma^ — MB)/Mp. The consis- 
tent picture that emerges for this splitting is sum- 
marized in a plot of a poster presented by Borrelli 
[37] at this conference. The mass dependence of 
this splitting is again rather weak, suggesting only 
an ~ 15% change between the b- and the c-mass. 
For the promising investigation of a number of 
other splittings, I refer the reader to [36]. 

7. CONCLUSIONS 

A significant advance has been made in the 
static approximation by applying the variational 
technique [4,5] to obtain ground state properties. 
Not only do refs. [4,5] obtain reliable numbers for 
the decay constants, but with the help of the gap 
computed in ref. [4] one can quantitatively esti- 
mate the contamination due to the first excited 
state that is present in other calculations. Soon, 
the precision of ^**''* may be limited mainly by 
the unknown 2-loop effects in the renormaliza- 
tion. 

Concerning the computations at finite mass, we 
point out that it is not sufficient to correct (ap- 
proximately) for one type of 0(o) effect or an- 
other. One needs to perform calculations over 
a range of lattice spacings with one action, one 
normalization of the fields, one definition of the 
meson mass and extrapolate to the continuum. 
The action of ref. [7] and/or the SW action should 
help in that they may allow for a smoother con- 
tinuum extrapolation than the standard Wilson 
action. We can also learn more about O(omy) 
effects once the 1-loop calculations of ref. [7] are 
finished. 

Higher-precision calculations are under way. I 
hope that by the time of the next conference the 



final errors on fjj and /g can be cut by about a 
factor 2 or more, especially if the variance reduc- 
tion described in section 5 is applied. 

In this review, we combined data from differ- 
ent groups to perform continuum extrapolations 
of certain quantities. As the reference scale we 
used the string tension because it is known with 
reasonable precision for the relevant range of /3- 
values. A related quantity, ro, is known to be 
much better for this purpose [38]. Once it will 
have been computed, we will not need to worry 
about the residual systematic errors in the deter- 
mination of the string tension. 

In order to finally compute the B — B mixing 
amplitude, the Ab = 2 four-fermi operator needs 
to be renormalized. It remains a true challenge to 
perform this renormalization non-perturbatively 
or "at least" to two loops. 
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